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BACKGROUND OF THE INVENTION 



I, Field of the Invention 

10 TECHNICAL FIELD 

The present invention relates to CDMA (Code Division 
Multiple Access) cellular telephone and wireless data, 
communications with data rates up to multiple Tl (1.544 Mbps) and 
higher (>100 Mbps) , and to optical CDMA. Applications are 
15 mobile, point-to-point and satellite communication networks. More 
specifically the present invention relates to a new and novel 
means — for a now approach to the design of waveforms and filters 
using mathematical formulations which generalize the Wavelet 
concept to communications and radar. 
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II, Description of the Related Art 



Multi-resolution waveforms used for signaling and/or 
5 filters and which are addressed in this invention are defined to 
be waveforms of finite extent in time and frequency, with scale 
and shift properties of multi-resolution over the time-frequency 
(t-f) space. These waveforms can also be referred to as multi- 
scale waveforms, and include multi-rate filtering and the 
10 Wavelet as special cases. The emphasis will be on digital 
design and applications with the understanding that these multi- 
resolution waveforms are equally applicable to analog design and 
applications. 

15 Background art consists of the collection of waveform and 

filtering design techniques which can, be grouped into six broad 
categories. These categories are: CI) least squares (LS) design 
algorithms for filters and waveforms that design to 
specifications on their frequency response, C2) analytic 

20 filters and waveforms which are specified by a few free design 
parameters that can be sub-categorized into current applications 
and primarily theorectical studies, C3) combinations of CI and 
C2 for greater flexibility in meeting communications and radar 
performance goals, C4) special design techniques . to control the 

25 noise levels from intersymbol interference (ISI) and adjacent 
channel interference (ACI) in the presence of timing offsets for 
multiple channel applications, C5) Wavelet filter design using 
scaling functions (iterated filter banks) as the set of design 
coordinates or basis functions, C6) filter and waveform design 

30 techniques for non-linear channels and in particular for 
operation in the non-linear and saturation regions of a high 
power amplifier (HPA) such as a traveling wave tube (TWT) or a 
solid state amplifier, and C7) LS dynamic filters derived from 
discrete filtering and tracking algorithms that include . adaptive 

35 equalization for communications, adaptive antenna, filters, 
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Wiener filters, Kalman filters, and stochastic optimization 
filters. 

Category CI common examples of LS digital filter design are 
5 the eigenvalue algorith m in "A New Approach to Least-Squares FIR 
Filter Design and Applications Including Nyquist Filters" [ 7 ] and 
the Remez-Exchange algorithm in "A Computer Program for Designing 
Optimum FIR Linear Phase Filters' . - HH~= — The eigenvalue algorithm 
is a direct LS minimization and the Remez-Exchange can be 

10 reformulated as an equivalent LS gradient problem through proper 
choice of the cost function. Both LS algorithms use the FIR 
(finite impulse response) digital samples as the set of design 
coordinates. LS design metrics are the error residuals in 

meeting their passband and stopband ideal performance as shown 

15 in FIG. 3. Category C2 common examples of analytical waveforms 
and filters for system applications are the analog Chebyshev, 
Elliptic, Butterworth, and the digital raised-cosine, and square- 
root raised-cosine. For theorectical studies, common examples 
are polyphase mul.tirate filters, quadrature mirror filters 

20 (QMF) , and perfect reconstruction filters. Although these 
theorectical studies have yet to yield realizable useful filters 
for system applications, their importance for this invention lies 
in their identification and application of ideal performance 
metrics for filter designs. Category C3, common digital example 

25 is to start with the derivation of a Remez-Exchange finite 
impuloo response — (-FIR-)- filter and then up-sample and filter with 
another bandwidth limiting filter. This results in an FIR over 
the desired frequency band that is larger than available with the 
Remez-Exchange algorithm and with sidelobes that now drop off 

30 with frequency compared to the flat sidelobes of the original 
Remez-Exchange FIR. A category C.4 common example is to select 
the free parameters of the category C3 filter in order to 
minimize the signal to noise power ratio of the data symbol (SNR) 
losses from the ISI, ACI, and the non-ideal demodulation. A r 

35 second common example is start with a truncated pulse whose 
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length is short enough compared to the symbol repetition 
interval, to accommodate the timing offsets without significant 
impact on the ISI and ACI SNR losses. This shortened pulse can 
then be shaped in the frequency domain. 

5 

Category C5 Wavelet filter design techniques discussed in 
the next section will serve as a useful reference in the 
disclosure of this invention. Category C6 common example is the 
Gaussian minimum shift keying (GMSK) waveform. This is a 

10 constant amplitude phase encoded Gaussian waveform which has no 
sidelobe re-growth through a non-linear or saturating HPA. 
Category C7 common examples are the adaptive equalization 
filter for communication channels, the adaptive antenna filter, 
and the Kalman filter for applications including target tracking 

15 and prediction as well as for equalization and adaptive antennas. 

Minimizing excess bandwith in the waveform and filter 
design is a key goal in the application of the C1,...,C6 design 
techniques for communications and radar. Excess bandwith is 

20 identified as the symbol a in the bandwidth- time product 
BT s =l+a where the two-sided available frequency band is B and the 
symbol repetition interval is T s . Current performance capability 
is represented by the use of the square-root raised cosine (sq-rt 
rc) waveform with a= 0.22 to 0.4 as shown in FIG. 48. The goal 

25 is to design a waveform with a=0 within the performance 
constraints of ISI, ACI, passband, . sideband, and passband 
ripple. This goal of eliminating the excess bandwidth corresponds 
to the symbol rate equal to the available frequency band 1/T S = 
B/(l+a) = B for a=0. This symbol rate 1/T S =B is well known to be 

30 the maximum possible rate for which orthogonality between 
symbols is maintained. A fundamental performance characteristic 
of ewthe new waveform designs is the ability to eliminate excess 
bandwidth for many applications. 
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Scope of this invention will include all of the waveform 
and filter categories with the exception of category C4 special 
design techniques and category C7 LS dynamic filters. Emphasis 
will be on the category C5 Wavelets to establish the background 

5 art since Wavelets are multi-resolution waveforms that can 
eliminate the excess bandwidth and have known design algorithms 
for FIR waveforms and filters. However, they do not have a design 
mechanism that allows direct control of the ISI, ACI, passband, 
sideband, and passband ripple. Category C2 theorectical studies 

10 also eliminate the excess bandwidth. However, they are not 

multi-resolution waveforms and do not have . realizable FIR design 
algorithms . So the emphasis in background art will be on 
Wavelets whose relevant properties we briefly review. 

Wavelet background art relevant to this invention consists 
of the discrete Waveform equations and basic properties, 
application of Wavelets to cover a discrete digital time- 
frequency (t-f) signal space, and the design of Wavelets using 
the iterated filter construction. Wavelets are waveforms of 
finite extent in time (t) and frequency (f) over the t-f space, 
with multi-resolution, scaling, and translation properties. 
Wavelets over the analog and digital t-f spaces respectively are 
defined by equations (1) and (2) as per Daubechies' s "Ten 
Lectures on Wavelets", Philadelphia: SIAM, 1992 
rof cronco — [1] . 

Continuous wavelet 

^,»('H«r'V(— ) (i) 

a 

Discrete wavelet 

(2) 

where the two index parameters "a,b" are the Wavelet dilation and 
translation respectively or equivalently are the scale and shift. 
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The \\f is the "mother" wavelet and is a real and symmetric 
localized function in the t-f space used to generate the doubly 

i t-l/2 

indexed Wavelet ij/ ab . The scale factor "|a| " has been chosen to 

keep the norm of the Wavelet invariant under the parameter change 
5 "a,b". Norm, is the square root of the energy of the Wavelet 
response. The Wavelets \|/ a , b and \\f are localized functions in the 
t-f space which means that both their time and frequency lengths 
are bounded. The discrete Wavelet has the time w t" replaced by 
the equivalent digital sample number "n" assuming the waveform is 
10 uniformly sampled at W T" second intervals. 

Wavelets in digital t-f space have an orthogonal basis that 
is obtained by restricting the choice of the parameters "a,b" to 
the values a=2~ p , b=qM2 p where *p,q" are the new scale and 
15 translation parameters and W M" is the spacing or repetition 
. interval T S =MT of the Wavelets (which from a communications 
viewpoint are symbols) at the same scale M p". Wavelets at M p,q" 
are related to the mother Wavelet by the equation 44+ 

20 

^ q (n) -2^^(2- p n-qM) (3) 



25 where the mother Wavelet is a real and even function of the 
sample coordinates. The orthonormality property means that these 
Wavelets satisfy the orthogonality equation with a correlation 
value equal to "1". 

XVp^m =1 iff both p = kandq = m 

n 

30 =0 otherwise (4) 
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Wavelet representation of a digital t-f space starts with 
selecting an N sample time window of a uniform stream of 
digital samples at the rate of 1/T Hz (1/second) equivalent to a 
"T" second sampling interval. The N point or sample t-f space 
5 in FIG. 1 illustrates a Wavelet representation or "tiling" with 
Wavelets that are designed analytically or by an iterated filter 
construction. 

The t-f space in FIG. 1 is partitioned or covered or tiled 

10 by a set of Wavelet subspaces tyV p ,p- O,l,—»0i-lj .1 where N=2 m . Each 

Wavelet subspace W p 1 at scale "p" consists of the set of Wavelet 
time translations {q = 0,1,..., N/2 P+1 -1} over this subspace. These 
Wavelet subspaces are mutually orthogonal and the Wavelets within 
each subspace are mutually orthogonal with respect to the time 
15 translates. This N-point t-f space extends over the time 
interval from 0 to (N-l)T 2 where — 3? — i« — tke — digital — oampling 
interval/ — and over the frequency interval from 0 to (N-l) 3 in 
units of the normalized frequency fNT 4. 

20 The iterated filter bank in FIG. 2 is used to generate the 

Wavelets which cover the t-f space in FIG. 1. Each filter stage 
5 consists of a high pass filter (HPF) and a low pass filter 
(LPF) . Output 6 of the LPF is . subsampled by 2 which is 
equivalent to decimation by 2. This t-f space space is an N- 

25 dimensional complex vector metric space V 7. At stage m in the 
iterated filter bank, the remaining t-f space V m _ x 8 is 

partitioned into V m+1 9 and the Wavelet subspace W m+ i 10 . 

Scaling functions and Wavelets at each stage of this filter 
30 bank satisfy the following equations 
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( Z»(n)=2- ,/2 Eh q ^(2n-q) Vq 

(5) 

Kn)=2 ,/2 Ig q9 )(2n-q) Vq 
q 

where 9 is the scaling function, \|/ is the Wavelet, HPF P 

coefficients are { h 9 Vq } , LPF P coefficients are { g g , Vq } , and 

5 the equations apply to the stages 0,l,...,m-l. Identifying the 
scale parameter and using the previous Wavelet formulations, 
enable these equations to be rewritten for stages p=0, 1, m-1 
as 

10 < 

q 

For eu ^the application the HPF P and LPF P are quadrature mirror 
filters (QMF) with perfect reconstruction. This means they cover 
the subspace V p with flat responses over the subband frequency 
15 including the edges of the frequency subband, and the HPF P 
coefficients are the frequency translated coefficients for the 
LPF P : {g q = (.-!)%, Vq}. 

Wavelet design using iterated filter bank starts with the 
20 selection of the scaling functions. Starting with a primative 
scaling function such as the one proposed by Daubechies — [1] y 
2_one can use the iterated filter construction given by equations 
(5) and (6) to derive successive approximations to a desired 
scaling function which has properties that have been designed 

25 into it by the selection of the filter coefficients {g^Vq} at 

each level of iteration. The Wavelets can be derived from these 
scaling functions using the iterated filter construction or 
scaling equations (5) and "(6) . 



(6) 
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Another use of the iterated filter construction is to 
design the scaling functions as Wavelets thereupon ending up with 
a larger set of Wavelets for multi-resolution . analysis and 
synthesis as illustrated by Coifman's Wavelets in "Wavelet 
5 analysis and signal processing" . - {-2+-= — Dcoign oxamploo abound in 
^he — open literature and the — recent publications — — HH-7 — [5] > 
-HH — arc roprcocntativc illustrations of tho diverge applications? 
of Wavelets. 

10 

SUMMARY OF INVENTION 

SUMMARY OF THE INVENTION 

15 This invention io a new and novel dcoign of Wavelet multi - 

re solution wavcf ormo and filters in the frequency domain with a 
p- eopcrty which provides a single waveform dcoign for all of the 
wavoformo at multiple rcoolutiono and ocalco. Wavclcto arc 
de oigncd to meet the opecific application rcquiromontO/ ar e 

20 complex, are defined by frequency harmonico, include a frequency 
^anolation parameter in addition to the ocalc and translation 
p aramctcro in current Wavclcto, and have the property that the 
ee t of dcoign coefficicnto in the frequency domain remain 
in variant for all poooiblc multi - resolution ocalcoy and time - 

25 frequency tranolationo . Frequency dcoign harmonico and the 
-f-E- cqucncy tranolation capability enable the waveform at multipl e 
r eoolution ocalco to be derived from the oinglo waveform dcoign 
by ocaling the dilation, tranolation, and frequency tranolation 
paramctcro uoing the oot of frequency dcoign harmonico. Wavel et 

30 waveform dcoign io illuotrated by Matlab 5.0 code to design — a 
linear filter waveform uoing an iterative leaot - oquarco 
ei genvalue approach to minimize the non - linear leaot oquaroo coo fe 
•g nnctionj and to ocalc thio waveform deoign for multi - rcoolution 
application opecified by the dilation, time tranolation, aftd 
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-£g cquoncy tranolation paramctcro. Additional rcoulto arc given 
£o r a conotant amplitude minimum - ohift - keying bandwidth-of f iciont - 
Hte dulation waveform and for a oynthotic aperture radar waveform. 

5 This invention ie — a — ftew — approach — fee — for the design of 

multi-resolution Wavelet w aveforms that improves their 

performance for engineering and scientific applications, and in 
particular for applications to communications and radar. Current 
practice is to 1) design traditional multi-resolution waveforms 

10 in the time domain using metrics which specify frequency 
performance and without consideration of Wavelet propertiesT — aftd 
or to 2) design Wavelet multi-resolution waveforms using the 
time-domain iterated multi-resolution filtering approach with a 
set of scaling functions used to perform the filtering and 

15 without direct considerations of the frequency performance. 
These two separate design approaches yield fundamentally 
different waveforms. Our now This invention provides a means to 
combine these two approaches to generate a new multi-resolution 
Wavelet w aveform with the best properties of the traditional 

20 multi-resolution waveforms and the Wavelet multi-resolution 
waveform. 

Omg ^This invention introduces %he — new— innovations for the 
design of etw?the now W avelet w aveform that 1) include the ISI 

25 (intersymbol interference) , ACI (adjacent channel interference) , 
and QMF (quadrature mirror filter) requirements into the design 
algorithm along with non-linear modifications to the traditional 
passband and stopband frequency requirements which are used in 
the Remez -Exchange and eigenvalue LS algorithms— -ana — rof oroncoo 

30 [8] , — R-] — respectively , 2) select the frequency harmonics as the 
design coordinates, 3) use an appropriate subset of the 
available Fourier domain frequency harmonics as the design 
coordinates with the property that this subset of coordinates is 
a basis for the multi-resolution Wavelet waveform design, and 4) 
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design the multi-resolution Wavelet waveform for no excess 
bandwidth a=0 . 

These innovations support the development of e ^the new 
5 waveforms 1) that are generalizations of Wavelets in the 
frequency domain using a means which makes them useful for 
communications and radar over their t-f space, 2) that can be 
designed to be an orthonormal basis or orthonormal set of 
coordinates over the available frequency bandwidth with the 

10 implicit property that the excess bandwidth vanishes a=0- to 
within the accuracies allowed by communication and radar design 
implementations, and 3) that have multi-scale properties which 
allow a single dc waveform design to be used to uniquely define 
the complete set of new waveforms over the t-f space for multi- 

15 scale applications Specific design algorithm examples developed 
in this invention disclosure are the LS design algorithms that 
use eigenvalue and gradient search techniques in the frequency 
domain to find the best waveform design which minimizes' the 
corresponding LS error residual cost function that is a weighted 

20 linear sum of the residual error metrics. 

Wavelet W aveform w aveform performance calculated in this 
invension disclosure from the application of these LS design 
algorithms, demonstrates the capability of this invention to 

25 provide a means to design waveforms that are improved over 
current practice. There are other algorithmic design algorithms 
which can be realized by modifications to this invention. An 
example given in this invention disclosure is the modification of 
these LS algorithms for application to constant amplitude 

30 bandwidth efficient (BEM) waveforms for communications. 
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BRIEF DESCRIPTION OF THE DRAWINGS AND 



THE PERFORMANCE DATA 

BRIEF DE S CRIPTION OF DRAWING S AND PERFORMANCE DATA 

5 The above-mentioned and other features,, objects, design 

algorithms, and performance advantages of the present invention 
will become more apparent from the detailed description set forth 
below when taken in conjunction with the drawings and performance 
data wherein like reference characters and numerals denote like 
10 elements, and in which: 

FIG, 1 is a Wavelet aft— N-point t-f space extending over the 
time interval (0, (N-l)T] and the frequency Interval (0, (N- 
. 1)/NT], which is tiled or covered by a set of orthonormal 
15 Wavelets at the scales p=0, 1, ...,m~l . 

FIG. 2 is a Wavelet iterated filter bank used to generate 
the set of Wavelets which tile or cover the t-f space in FIG. 1. 

20 FIG. 3 illustrates the power spectral density (PSD) of 

the new Wavelet waveform at dc, and the set of stopband and 
passband design requirements. 

FIG. 3 io tho power spectral dcnoity (PSD) — template for the 
25 tho stopband and paooband uocd to construct the LS error metrics, 

FIG. 4 is the flow diagram of the LS metrics and cost 
functions and the final cost function used to find the optimal LS 
30 solution for our the Wavelet waveform. 

FIG. 5 is the flow diagram for the hS — eigenvalue recursive 
solution algorithm to find the optimal harmonic coordinates that 
define the LS solution for emg -the Wavelet waveform. 
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FIG. — 5 io tho Matlab 5.0 code for the LS algorithm used to 
design tho Wavelet waveform in FIG. 6 . 



FIG. 6 plots the dc PSD in dB units for the frequency 
response of the mother Wavelet at dc and the square-root raised- 
cosine (sq-rt r-c) waveforms with excess bandwidth parameter a= 
10 0.22, 0/40 vs. the normalized frequency in units of symbol rate 
1/T 8 . 

FIG. 7 plots the dc PSD in dB units for our waveformo vs. 
the normalized frequency fT s for the f*ew — Wavelet waveform 
15 designed with modifications to the LS algorithms for application 
to constant amplitude BEM communications, and for a Gaussian 
minimum shift keying (GMSK) waveform. 

FIG. — 7 — plots — the — de — PS© — A** — dB units — So* — the — frequency 
20 response — of the mother Wavelet at dc and an, optimized Gaussian 

' minimum — shift — keying (GMSK) waveform =vs-= the — normalized 

frequency fTV =- 

FIG. 8 plots the amplitude of the dc radar ambiguity 
25 function vs.' the normalized frequency fT p and normalized time tT c 
where the pulse time T p and chip time T c are both equal to the 
communications symbol interval T s for this example, for the new 
waveform and an unweighted chirp waveform. 
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DISCLOSURE OF THE INVENTION 
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Wavclot multi - resolution waveforms incorporate a frequency 
tranolation, arc complex/ arc designed in frequency harmonic 
coordinates, a oinglc mother Wavelet design io used for all 
scales^ — and in the t - f apace arc realized ao FIR filters and 
5 waveforms . 

New Thcoe Wavelet waveforms in this invention disclosure are 
generalizations of Wavelets in t-f space which enable them to be 
useful for communications and radar applications* This 

10 generalization is accomplished by 1) the addition introduction 
of a frequency translation, by 2) relaxing — changing the 
orthonormality condition in equation (4) to apply to waveforms 
within the same space {q} and over the scales {p} with the 
inclusion of the frequency translation, and 3) by their 

15 characterization and design in the Fourier domain. With frequency 
translation the analytical formulation of these new waveforms as 
a function of the baseband or mother waveform centered at dc( dc 
refers to the origin f-0 of the frequency space) becomes 

20 

New waveform ao a function of dc waveform 

¥pqr (n) = 2~ pa K2~ p n - qM ) e i2/rfc (p ' r)nT ( 7 ) 

25 where f c (p,r) is the center frequency of the frequency translated 
dc waveform, at scale M p" and frequency index w r". The purpose 
of the frequency index "r" is to identify the center frequencies 
of the waveforms at the scale "p" in the t-f space. The dc or 
baooband or mother waveform is the generalization of the mother 

30 Wavelet for multi - resolution waveforms. 

These waveforms satisfy the complex orthonormality 
equations 
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Orthogonality equations 

Zv =1 iff p = kandq = m and r = v 

£ Y p,q,r Y k,m,v 

= 0 otherwise (8) 

where is complex conjugation, and are generalizations of the 

5 orthonormality equations for the analytical Wavelets in (4) . 

The new waveforms in equation (7) expand the Wavelet 
analytical formulation to include a frequency variable. 
Wavelets are functions of the scale and translation parameters 

10 xx p,q". The concept of an additional parameter to provide an added 
degree of flexibility in their tabulation was first introduced by 

the — Coif man — library — ef Wavelets — — which — tree — a** — oocillaton 

parameter £er tabulation — that roughly corresponds — fere the 

frequency of oscillation of the Wavelets. — The frequency variable 

15 together with the Fourier domain design are entirely new means 
for deriving these new waveforms as generalization of the 
traditional Wavelets aftd — their — modification — with — the — added 
. parameter in — [2] . 

20 New —These waveforms are generalizations of Wavelets in the 

.frequency domain using a means which makes them useful for 
communications and radar in the t-f space. The basis vectors for 
this metric space V consist of a subset of the admissible set 
of scaled and- translated waveforms. {^ Mr ,Vp,q,r} derived from the 

25 dc waveform y/ as per equation (7) . An admissible waveform is 
any combination that covers F=t-f space. We are interested in 
the Fourier domain representation of the dc waveform \f/ in V , 
and in particular in a subset of the discrete Fourier transform 
(DFT) harmonic coefficients over the Fourier domain which we 

30 intend to use as the design coordinates. Starting with the z- 
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transform and continuous Fourier transform, the DFT harmonic 
coefficients. are defined by the following equations 



DFT Harmonic Coefficients (9) 
W( z ) = £<K n ) z*° z- transform 

n 

y/(co) = Z^(n) e' l<wl Fourier transform 

n 

Vv. = V<& 

= 2>(n) W^" DFT harmonic coefficients for Vk 

where 

iff (co) = (1/ N») J> k £ e i(2 ^ ^ 

k n 

= £yV sin( (*>/2 - ak/N 1 )^) / Nsin(<a/2 - ;rk/N') 
k 

= [Harmonic interpolation for" k"] 

k 

{k} = DFT frequency or harmonic coefficients such that 4 = k/N f T 

where f k is the harmonic frequency corresponding to "k" 
N f = length of ^(n) 

where ^(a)) = \y/(a))f is the power spectral density of the dc 

waveform. These equations define the frequency representation of 
o** gthe new waveforms in terms of the available set of harmonic 
coefficients, which set is considerable larger than required for 
15 most applications. 

** ew — The new waveforms are an brthonormal basis with no 
excess bandwidth, — et&e — propcrtico — that — we — will — demonstrate . 
Orthomality and no — oxcooo bandwidth arc which p roperties which 
20 are asymptoticly approached by ea ^the new waveforms to within 
design accuracies inherent in communications and radar. To 
proccod — demonstrate these properties w e need to identify the 
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structure of the dc waveform in F. We start with definitions 
for the parameters and coordinates in the following equation 
(10) . The waveforms derived from the dc waveform- will be 
designed to be orthogonal over both time translates "MT" and 
5 frequency translates "1/£^MT" which respectively correspond to 
the Wavelet symbol spacing T S =MT and the adjacent channel 
spacing 1/T S . This means the orthogonal spacing of the waveforms 
in V are at the time-frequency increments (MT, 1/irfMT) = (T s . 
1/T S ) . In the interests of constructing e ^the aew— orthonormal 

10 multi -resolution waveforms to cover V it will be convenient to 
assume that M,L are powers of 2. We need the following 
definitions for the parameters and coordinates. 

15 



20 



25 



30 
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Parameters and Coordinates 



(10) 



N r = Length of xj/ which is an even function about the center and 
which spans an odd number of points or samples 
= ML + 1 where M, L are assumed to be even functions for 

convenience of this analysis 
= Number of points of v 
M = Sampling interval for y/ 

= Spacing of y/ for orthogonality 
L = Length of \y in units of the sample interval M 

= Stretching of y over L sample intervals 
n = n 0 + njM 

= partitioning into an index n 0 over the sample 

lengthn 0 = 0,1,.. .,M- land an index n, = 0,1...,L-1 
over the sample intervals 
k = k 0 + k,L 

= partitioning into an index k 0 over the harmonic 
frequencies k 0 = 0,1 ,...,L - 1 corresponding to the 
stretching and an index k x = 0, 1 , . . . , M over the harmonics 
frequencies corresponding to the admissible frequency 
slots for y/ 

The harmonic design . coordinates are selected using the 
following observation. For most applications and in the following 
development it is assumed that the waveforms are spectrally 
^contained in the frequency interval l/3x£ — MT corresponding to the 
frequency spacing. This suggests the harmonic design coordinates 
be restricted to the. subset of L harmonics {k 0 =0 f 1, . . . ,L-1} 
covering this spacing. These L harmonics correspond to the 
stretching of the mother waveform over the L repetition 
intervals. 

Obviously, for some applications as will be demonstrated 
later, the spectral containment is spread out over several 1/LT 
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frequency increments whereupon one must increase the subset of 
design harmonics to possibly 2L, 3L or larger. 



Pursuing this application uoing L dooign harmonics^ — The the 
5 DFT equations for the dc waveform in (9) when rewritten in terms 

of the L harmonic design coordinates { y/ k , Vk 0 } become: 
DFT equations for dc waveform 

10 

= /^V( n ) W^° n harmonic design coordinates 

n 

y/(n) =(1/N')5;^ k W*° n new waveform defined in terms (11) 
of the L harmonic design coordinates {^ o ,Vk 0 } 

15 

It will now be shown that It will bo proven that the use of 
these L harmonic design coordinates is sufficient -fe* — e*H? — F*ew 
wavof orm — under time translates to be a basis for the 
corresponding subspace of V— This w hich m eans these waveforms 
20 provide a complete set of coordinates to describe this subspace. 
We will uoq the 

will use the theorems of Karhunen-Loeve and Mercer and will 
limit the demonstration to the dc waveform for simplicity and 
without loss of generality. We start by considering Consider the 

25 expansion of a random complex sequence {z(n), Vn} in a series of 
waveform coordinates consisting of time translates of the mother 
waveform y/ . The sequence {z(n) , Vn} is a zero-mean stationary 
random process which is orthonormal over the sample interval "M" 
and has a frequency spectrum which is flat and extends over the 

30 frequency range 1/LT which is centered at baseband corresponding 
to a zero frequency. This means the { z (n) , Vn} cover the 
subspace of V corresponding to the scale of e ^the ftew — dc 
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waveform y/ and its time translates { y/(n-qM) = if/ q (n\ Vq } . In 

addition, the V is now considered to be extended over a time 
interval which is relatively large compared to the N-dimensional 
t-f space in FIG. —1,2 to avoid end-effects on the analysis. We 
5 start by approximating the sequence {z(n), Vn} by the { z(ri) , Vn} 
where : 

z(n)=2:Z q Kn-qM) 

= £Z q ifin 0 + (n,-q)M) 

q 

where the complex coefficients {z q }are 

derived from the original sequence using the ( 12 ) 

orthogonality properties of y/ 
Z q =Xz(n)^*(n-qM) 

n 

= Iz(n)^(n) 

n 

10 

The following equations prove that the coefficients {Z q , Vq} are 
orthonormal : 

Z q Z q .* = Ez(An-qM)z*(An-q , M)^ q (An)^ q . *(An) (13) 

An 

= ^qq« E|^ q (n)| 2 since the sequence {z} is orthonormal 

n 

15 for qM time translates 

= S , with normalization of the energy of if/ 



Equations (12) and (13) together prove the Karhunen- 
Loeve's theorem which proves the following equation for the 
20 accuracy in approximating the stochastic sequence {z(n), Vn} by 
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Vn}. This accuracy is expressed by the expected xx E(o)" 
squared error w (o)" in this approximation: 

\ ■ ■ 

5 E{z(n) - z(nf } = 1 - 5>(* + n x M)\ 2 (14) 

We need to prove that the right hand side of this equation is 
zero which then proves that the approximating sequence is equal 
to the original sequence in the mean-square sense. In turn this 
10 proves that the new waveform coordinates { y/ q , Vq } are a basis for 

the original sequence {z(n), Vn} which is o^ ^the goal. 

The right hand side of equation (14) when set equal to zero 
expresses Mercer's theorem so etw?the goal is to prove Mercer's 
15 theorem. To do this we use the DFT of y/ in equation (11) and the 
coordinates in (10) to evaluate the right' hand side of equation 
(14) « We find 

where T = (l/Lj^Wp* 

= sin(;r Ak 0 )/Zsin(;r Ak 0 1 L) 

= 1 for Ak 0 = 0 (15) 

= 0 otherwise 
This proves that 

1- +n x M)\ 2 = 0 Vn 0 

20 

which proves that equation (1&14) reduces to 
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E{z{ri)-z{nf\= 0 



which as per the above proves that the set of multi-resolution 
5 waveforms is a basis. This proof easily generalizes to the 
multi-resolution waveforms at all of the scales {p} and time 
translates {q} , and to the expansion of the harmonic design 
coordinates over 2L, 3L, ... as required by the application. 

10 New — Next we demonstrate that these waveforms have multi- 
scale properties^ that we will demonstrate. We start , with the 
observation that the Tfee — Fourier domain design for the new 
waveforms provides a natural and easy way to derive the complete 
set of waveforms {^ ?5 Vp,q} for the space V from the design of 

15 the dc waveform iff by using the same invariant set of Fourier 

domain harmonic design coordinates {-^ 9 Vk 0 } derived for the dc 

multi-resolution waveform. This demonstration requires that we 
show 1) how the multi-scale transformations are implemented with 
the design in the Fourier domain, and 2) how the waveform design 
20 remains invariant, under scale changes. 

First consider the multi-scale transformation which derives 
the waveforms at the scale and shift parameters "p,q" for from 
the dc waveform at scale "p=0" and centered at the origin x> q=0". 

25 We begin by extending the parameters and coordinates in equation 
(10), to include both scaling and . subsampling or decimation in a 
form that is equivalent to the. iterated filter bank construction 
which is used to derive the current Wavelet waveform using the 
filter scaling functions. Starting with the coordinates at scale 

30 "p=0" the parameters and coordinates at scale "p=p" are given by 
the equations: 
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5 ^E = 0 parameters and coordinates 

n = n 0 + rii M (17) 
n 0 = a 0 + a x 2 + ... + a m -i 2 m_1 ' 

= M points 
M = 2 m . 

10 n x = b 0 + b x 2 + ... +bi_i 2 1 " 1 

= L points spaced at M sample intervals 
L - 2 

p=p parameters and coordinates 

15 2* p n(l2 p ) = scaled by w 2" p// 

and subsampled or decimated by 2 P :1 
= n 0 (p.) + ni(p)M 

2" p n 0 (l2 p ) = scaled by "2" p " 

and subsampled or decimated by 2 P :1 
20 = n 0 (p) 

= ap + ap_i2 + ... + ap+m-i 2 P 

= M points spaced at 2 P sample intervals 

. 2" p n x ( I 2 P ) = scaled by u 2' p " 

and subsampled or decimated by 2 P :1 
25 = m(p) 

= b p + b p+ i 2 + ... + b p+1 -! 2 P+1 " X 

= L points spaced at M2 P sample intervals 



30 together with the observation that the sampling interval "T" is 
increased to W 2 P .T" under the scale change from "p=0" to "p=p" and 
subsampling or decimation from "1:1" to "2 P :1". Combining these 
equations with the analytical formulation in (7) and the Fourier 



23 



domain representation in (11) enables the waveforms at the 
parameters "p,q" to be written as a function of the Fourier 
domain harmonic design coordinates: 

5 

¥pjtf («) = 2- pQ ^(2" p n - qM ) e arf < (p ' r)nT 

= (2 -W2 /N<) £^ w m^m, (18) 

for all admissible scale, translation, and frequency index 
10 parameters "p,q,r". 

Next we need to demonstrate that the frequency domain 
design in (11) remains invariant for all parameter changes and in 
particular for all scale changes.. This multi-scale property 

15 expresses the accordion behavior of the design in that the 
Wavelets at . different scales are simply the stretched and 
compressed versions of the mother waveform with the appropriate 
frequency translation indices. This multi-scale invariancy means 
that the design for a M=16 channel filter bank remains the same 

20 for M=100 or M=10,000 channel filter banks, when the overlap L 
and the performance goals remain constant. To demonstrate this 
invariant property across scales, we consider the hh? — multi- 
resolution waveform at scale xx p" with the other parameters set 
equal to zero for convenience >x q-0, r=0" and without loss of 

25 generality. The Fourier domain frequency response *F(/) can be 
evaluated starting with the original formulation in equation (7) : 



30 
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DFT at "p, q=0, r=0" (19) 



%(/) = (i/N')2X Y,w-w' 2 ' T *< 

k 0 n(p) 



)n(p) 



k 0 



sin(;r (fN'2 p T-k 0 )) 



N'sin(;r (fN7 p T-k 0 )/N')^ 
= [Harmonic interpolation for "k 0 M ] 



10 This only differs from the harmonic representation in equation 
(9) in the restriction of the design coordinates to the subset of 

harmonic coefficients { y/ k , Vk 0 } and the stretching of the time 

interval to «2 P T" corresponding to the scale "p=p". The harmonic 
interpolation functions are observed to remain invariant over 

.15 scale changes upon observing that the * frequency scales as 
"f-I/2 p T" which means the w f requency*time" product remains 
invariant with scale changes as per the fundamental property of 
the waveforms. This means the harmonic interpolation functions 
remain invariant with scale change and therefore the frequency 

20 response remains an invariant. This demonstrates the waveform 
design is an invariant across the waveform scales which means we 
only need a single design for all scales or resolutions of 
interest, 

25 LS design algorithms for new waveform will be described to 

illustrate the advantages eua ethe new w aveform has over current 
designs. The two LS algorithms described are the eigenvalue and 
the gradient search which respectively can be reduced to 
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algorithms which are equivalent to J fehe original — current 

eigenvalue -FH — and Remez-exchange -HH — waveform design algorithms 
for application to a uniform filter bank. We consider the t-f 
space which is spanned by a uniform polyphase filter bank 
5 consisting of M channels at the frequency spacing f w =l/MT where T 
is the digital sampling interval, and the filter waveform FIR 
time response is stretched over L sampling time intervals T s . 
This polyphase filter bank is ideally decimated which means the 
filter ouput sample rate 1/T S is equal to the channel-to-channel 

10 spacing T S =MT, equivalent to stating that there is no excess 
bandwidth a=0. O ^The design for this topology is immediately 
applicable to an . arbitrary set of multi-resolution filters 
through the scaling equation (18) which gives the design of 
omg -the waveform at arbitrary scales in terms of e ^the design of 

15 the dc waveform. 

For this polyphase filter bank used to construct the dc 
waveform or filter impulse response, e ^the LS example design 
algorithms will use 5 metrics consisting of the 2 prior art 

20 passband and stopand metrics, and the 3 new metrics consisting 
of the ISI, ACI, and QMF, and solve the LS minimization problem 
using as design coordinates the subset of harmonic coordinates 
which are a basis. Since ewthe 2 — two example LS design 
algorithms only differ in the use of an eigenvalue LS 

25 optimization and the use of a gradient search LS optimization, 
the flow diagrams for the construction of the cost functions and 
the solution for the optimal waveform will . be identical. 
However, there are differences in the construction of the cost 
functions from the respective metrics and in the iterative 

30 solution mathematics. Both LS solutions for the harmonic 
design coordinates minimize the weighted sum of the error 
residuals or cost functions from the 5 metrics. These design 

coordinates — are the Fourier harmonics {^,Vk 0 } for {k 0 = 

0, 1,...,L-1} . Resulting algorithms are easily extended to the 
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applications requiring the design coordinates to cover 2L, 3L, ... 
harmonics • 

Frequency domain design coordinates are related to the 
5 waveform time domain digital samples or coordinates as follows. 

Mappings of timo — ^ — frequency to time . (20) 

Time domain design coordinates {y(w),Vn}are real and 

symmetric and can be represented by the reduced set {h t (n), 
n=0,l,...,ML/2} 

h t (n) = j/(0) for n=0 

= 2y/(n) for n=l, 2, ...,ML/2 
= time domain design coordinates 

Frequency domain harmonic design coordinates {^ 5 Vk 0 } are 

real and symmetric and can be represented by the reduced 
set {hf(k), k=0,l, ,...,L-1} 

h f (k) = yr k for k=k 0 =0 • < 

= 2y/ ko for k^ko^l^.^L-l 
= frequency domain design coordinates 

Mapping' of the frequency coordinates {h f (k), k=0, 1,...,L-1} 
into the time coordinates {h t (n), n=0, 1, ...,ML/2} is defined 
25 by the matrix transformation to within a scale factor 

h t = B h f 

where 

h f = (hf(0) r ...,h f (ML/2) ) t transpose of column vector 
h t = (h t (0), ...,h t (ML/2) ) t transpose of column vector 
30 B = (ML/2 +1) x L matrix 

= t Bkn ] matrix of row k and column n elements B kn 
Bkn = 1 /ML for n=l 

— '= 2 cos(27ikn/ ML) otherwise 

27 
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Mappings of the — frequency coordinates — R Mk) > — k-Oy 1,...,L - 1} 
into — ^febe — time — coorciinatoo — H*t (n) / — n-0, 1, ...,ML/2] — attd — vice 
vcroa arc defined by the matrix transformations 
k^_=-B — h s 



where 

h*-= — fh *(0) > ...,h^(ML/2) ) — tranopooo of column vector 

10 : = — H MO) / ...,h % (ML/2) )* — tranopooo of column vector 

B — ~ (ML/2 + 1) x L matrix 

— : ~ [ B^H — matrix of row k and column n olomento EW 



B te "" 1 /ML 1 1 for n-1 

- 2 coo(2nkn/ ML + 1) othcrwioc 

15 B7 * - L x (ML/2 1 1) matrix 

= — HB^S matrix of row k and column n elements - B-g*- 

-B^=— 1 ■ for n-1 

: = 2 coo(2rikn/ ML + 1) othcrwioc 



and 



20 B~*~B - I L x L identify matrix 

: — B— B - ^-^— I (ML/2 i l) x (ML/2 H) identity matrix 

wherein the N /= =ML+1 has been replaced by ML since a single end 
point has been added to the the FIR to make it symmetrical for 
25 ease of implementation for the example Wavelet being considered 
with a sample at the mid-point that makes the number of samples 
N' an odd number. 



30 Passband and stopband metrics and cost functions are 

. derived with the aid of FIG. 3 which defines the po wer opcctral 
density — (-PSD-^ parameters , of interest for the passband and 
stopband of the PSD *F(6>) for communications applications. 
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Requirements for radar applications include these listed for 
communications* Referring to FIG. 3 the passband 11 of the 
waveform PSD is .centered at dc (f=0) since we are designing the 
dc or baseband waveform, and extends over the frequency range C0p 
5 extending from -G>p/2 to +0)p/2 12 in units of the radian 

frequency variable (0=27ifT 13 where T io the digital sampling 

interval defined in FIG. 1 . The frequency space extends over the 
range of f=-l/2T to f=+l/2T which is the frequency range in 
FIG 1 translated by -1/2T so that the dc waveform is at the 

10 center of the frequency band. Quality of the PSD over the 
passband is expressed by the passband ripple 14. Stopband 15 
starts at the edge 16 of the passbands of the adjacent channels 
+ /-o> a /2 16 and extends to the edge of the frequency band ©=+/- 
n 17 respectively. Stopband attenuation 18 at +/-co a /2 

15 measures the PSD isolation between the edge of. the passband for 
the dc waveform and the start of the passband for the adjacent 
channels centered at +/-o) s 19. Rolloff 20 of the stopband is 
required to mitigate the spillover of the channels other than the 
adjacent channels, onto the dc channel. Deadband or transition 

20 band 21 is the interval between the passbands of contiguous 
channels, and is illustrated in FIG. 3 by the interval from <Dp/2 
to (o a /2 between the dc channel and adjacent channel at co a . 
Waveform sample rate G) s 22 is the waveform repetition rate. 
For the LS example algorithms, the waveform sample rate is equal 

25 to. the channel-to-channel spacing for zero excess bandwidth. 
Therefore, 1/T S = <o s /2rcT=l/MT which can be solved to give a> s =27i/M 
for the radian frequency sampling rate of the filter bank which 
is identical to the waveform repetition rate. 

30 We start by rewriting the DFT equations for the dc waveform 

in (11) as a function of the {h f (k), k=0, 1,...,L-1} 
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y/{o>) = ]^^(n)cos(nfi>) 

n 

= c* B h f using (21) and the definition of the vector "c" (21) 
c = ( 1, cos(#>), cos( (ML/2) a>) ) % transpose of column vector 

which is equivalent to the equation for y/{(o) in (9) expressed as 
5 a linear function of the {h f (k), k=0, 1,...,L-1} . However, — this 
functional — form ouggcotcd by — the — eigenvalue — formulation — R-] — As- 

more — convenient — to — analyze . An ideal "c" vector "c r " will be 

introduced for the passband and the stopband in FIG. 3-r in order 
to identify the error residual 8y/(oj) at the frequency "m" in 
10 meeting the ideal passband and stopband requirements. The ideal 
PSD is flat and equal to w l" for the passband, and' equal to "0" 
for the stopband. We find 

15 Error residuals for passband and stopband (22) 

c r = (1,1, .„ r l) t passband ideal w c" 
= {0 f 0 f ^,0) t stopband ideal M c" 
8c = c r - c error vector 

2Q 8 ¥ {a>) = & l Bh f 

= Residual error in meeting the ideal spectrum at "a>" 

The LS metric for the passband and stopband can now be 
constructed as follows for the eigenvalue and the LS optimization 
25 (or equivalently, the LS algorithm) design algorithms— 
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( 



Passband and stopband metrics 



(23) 



J (band) = — - — (|<^(<2>)| 2 do) Eigenvalue 
band J 



band 



= \Sy/ 



Eigenvalue 
LS 



(23) — continued 



-where 



10 



15 



band = [0, Op ) 
= '{<o s , n] 



passband 
stopband 



— |(B t && t B)d6) 



band band 

- L x L matrix 

8v|/ = ( 8v[/((Oi) 8\|/(Ou) ) t 

- vector of error residuals at the 
frequencies ©i , a>u across the band 

= norm or length of the vector (o) and which 

includes a cost function for the errors of 
the individual components 



20 



25 



where it is observed that the eigenvalue approach requires that 
the LS metrics be given as quadratic forms in the design 
coordinates {h f (k), k=0, 1, L-l } whereas with the LS approach it 
is sufficient to give the LS metrics as vector norms with 
imbedded cost functions or ah equivalent formulation. 
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QMF metrics express the requirements on the deadband that 
the PSD' s from the contiguous channels in FIG. 3 add to unity 

across the deadband [(Op , <o s ] in order that the filters be QMF 

filters. By suitable modification of the error vector 8c, the 
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previous construction of the passband and stopband metrics can be 
modified to apply to the deadband. This la a Wavelot requirement 
eft — the — frequency — coverage — of the — Wavelet basis — a-s — observed — if* 
Table 1. We find 



10 
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Deadband metrics 



J (deadband) = 



1 



deadband 



§Siff(ct))\ do 



deadband 



where 



8c = c r - c((D) - c(tc/M - ©) 



(24) 

Eigenvalue 

(24) continued 
Eigenvalue 

LS 



where c((D)=c as defined in (22) and (23), and c(n/M - o) = c at 
20 the offset frequency "n/M - co" corresponding to the overlap of 
the contiguous filters over the deadband. 



Orthonormality metrics measure how close we are able to 
designing the set of waveforms to be orthonormal over the t-f 

25 space, with the closeness given by the ISI and the ACI. ISI io 
%fee — non - orthogonality — error — between — channel — output — samples 

separated by multiples — &4 — the — sampling — interval 1/MT — seconds 

where T io the — sample time and M io the interval of contiguous 
samples . — AGI — is — the — non - orthogonality — error — between — between 

30 channel — output — samples — within — a — channel — a**d — the — samples , in 
adjacent — channels — at — the — same — sample — time — and at — sample — times 
separated by multiples of — the — sample — interval . Ao observed ao 
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noioc contributiono within each oamplo in a givon channel > fefee 

ISI io the noioc contribution duo to the other received waveforms 
at the different timing of f oeto correoponding to multiploo of the 

oampling interval. Likewise, tho ACI io tho noioc contribution 

5 €tee — — the — other — wavcf ormo — in adjacent — channelo — a* — fehe — oamc 
oampling time and at multiploo of the oampling interval. 
ISI and ACI errors are fundamentally caused by different 
mechanisms and therefore have separate metrics and weights to 
specify their relative importance to the overall sum of the LS 

10 metrics. ISI is a measure of the non-orthogonality between the 
stream of waveforms within a channel as per the construction in 
FIG. 3. On the other hand, ACI is a measure of the non- 
orthogonality between the waveform within a channel and the other 
waveforms in adjacent channels. This means the stopband 

15 performance metric has a significant impact on the ACI due to the 
sharp rolloff in frequency of the adjacent channel, and the ACI 
metric is then a measure of the residual non-orthogonality due to 
the inability of the stopband rolloff in frequency from 
completely eliminating the ACI errors. 

20 

We assume that the received waveform is identical to the 
filter waveforms and is transmitted at the filter output sample 
intervals equal to -t/MT seconds. The second assumption means we 
are assuming the receiver is synchronized with the received 
25 signal. Since there is no information lost by sampling 
asynchronously with the received waveform, we are free to make 
this synchronization assumption without loss of generality. ISI 
metrics are derived in the following set of equations. 

30 ISI metrics (25) 

Mapping of h f into \\? 

\|/ = (\|/(-ML/2), ... , \\f{ML/2)) t transpose of column vector 
= H h f 
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H = (ML+1) x L matrix of elements 

H kn =1 for n=l 

= 0.5 cos (27ikn/ML+l) otherwise 

5 Offset matrix A 

A = L x (ML+1) matrix of elements A kn 

Akn = 0 for row vector k=l 

- [0 0 ... V|/(-ML/2) ... V|/( (L-k)M+l) ] row, vector k 
10 (25) continuoci 

12 kM LM+1 column elements 

ISI error vector 6E 
8E = L x 1 column vector 
15 = A H hf 

IS I metric 

J(ISI) = 8E fc 8E Eigenvalue 

= Non-linear quadratic function of h f 

20 = \\SE\\Z ls 



ACI metrics are derived using the ISI metric equations with the 
following modifications. 
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ACI metrics (26) 



Mapping of h f into \\t is the same as developed for ISI 
Offset matrix A elements are changed as follows to apply to 
30 Channel 1: 

A^ - 0 for row vector k=l 

= [0 0 ... \|/(-ML/2)^ ... V ( (L-k)M+l) W£- k) ] 
1 2 kM+1 LM+1 
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which means the ACI error vector SE is 
8E = L x 1 column vector 

= A H he : h f 



5 ACI metric for the two contiguous channels 

J(ISI). « 2 SE t 8E Eigenvalue 

= Non-linear quadratic function of h f 

= 2 \\5E\\Z LS 

10 where the factor "2" takes into account there are two contiguous 
channels or one on either side of the reference channel 0 in FIG. 
3. Because of the fast rolloff of the frequency spectrum the 
addition of more channels into the ACI metric is not considered 
necessary, although the functional form of the ACI metric in (26) 

15 allows an obvious extension to any number of adjacent channels 
which could contribute to the ACI. 

Cost function J for the LS algorithms is the weighted sum 
of the LS metrics derived in (23) , (24) , (25) , (26) . The LS 
20 algorithms minimize J by selecting the optimal set of frequency 
coordinates {h f (k),Vk} for the selected set of parameters used 
to specify the characteristics of the dc waveform, frequency 
design coordinates, LS metrics, and weights. Cost function and 
optimization techniques are given by the equations 

25 

Cost function J (27) 
J = ^w(metrics) J(metrics) 

X . 

= weighted sum of the LS metrics J (metrics) 

30 where 

metrics = passband, stopband, deadband, ISI, ACI 
{w (metrics) , ^/metrics} - set of weights 
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y^w(metrics) = 1 normalization 

X 



Optimization goal 

Goal: minimize J with respect to the selection of 
the {h f (k),Vk} 

Optimization algorithms 

Two algorithms are the Eigenvalue and the LS optimization 
where the eigenvalue optimization algorithm uses the non- 
linear quadratic formulations of the LS metrics and the LS 
optimization algorithm uses the norm formulations for the 
LS metrics, 

FIG. 4 is a summary of the LS metrics and the construction 
of the cost function J for LS algorithms which use quadratic 
error metrics and LS algorithms which use norm-squared metrics . 

Design parameters 23 are the input and output design 

parameters. Input parameters are the number of polyphase 
channels M or equivalently the number of digital samples at 
spacing T over the symbol interval T S =MT, the length of the FIR 
time response for the waveform in units of L which are the number 
of' digital samples per waveform repetition interval T s so that 
the total number of digital samples for the symmetric FIR time 
response is equal to N^ML+1, number of DFT samples per FIR 
length n_fft for implementation of the LS algorithms, passband 
radian frequency G>p, stopband radiam frequency co a , waveform 
repetition rate in radian frequency G) s , selection of the set of 
design coordinates {h f } to be used in the optimization, and the 
metric weights {w (metrics) } . Output parameters are the set of 
harmonic design coordinates {h f } that minimize J. Band metrics 
24 are the passband, stopband, and deadband metrics defined in 
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equations (23), (23), (24) respectively. Interference metrics 
25 are the ISI and ACI metrics defined in equations (25) and 
(26) respectively, LS cost function J 26 is the weighted 
linear sum of the metrics defined for the band 24 and the 
interference 25 as defined in equation (27) 

FIG. 5 is a flow diagram of the LS recursive solution 
algorith m using quadratic error metrics and norm-squared error 
metrics. There are two loops with the topology constructed so 
that the outer loop 27 is an iteration over the set of metric 
weights {w (metrics) } , and the inner nested loop 28 is the 
recursive or iteratived LS solution to find the optimal {h f } for 
the given design parameters and weights where each step refers to 
an iteration step in the solution. A recursive LS solution is 
required to be able to solve the highly non-linear interference 
metrics as well as the band metrics when one chooses to use non- 
linear techniques to construct these metrics. 

Consider the inner nested loop 28. The recursive LS 
solution starts with the initial step i=0 29 which begins with 
the selection of the input design parameters and weights and the 
selection of an initial set of values for the {h f } 30. Next the 
band metrics are calculated 31. Since this is the initial step 

i=0 32 the cost function J is restricted to the linear band 

V. 

metrics in order to fine the linear approximation to {h f } 33 to 
initialize the non-linear solution starting with step i=l 34 
wherein the highly non-linear interference metrics are included 
in J. Following the signal flow, for step i=l the band metrics 
31 and interference metrics 35 are calculated, weighted, and 
summed to form the cost function J 36. An iterative solution 
algorithm 37 finds the best approximation in step i=l to the 
{h f } which minimizes J using the approximation to {h f } from the 
previous step i-0 to linearize the search algorithm coefficients 
for the eigenvalue and LS step i-1 iteration. If there is no 
convergence to the correct solution for {h f } this recursive 
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for the eigenvalue and LS step i=l iteration. If there is no 
convergence to. the correct solution for {h f } this recursive 
sequence of calculations is repeated for step i=2 34. This 
recursive solution technique is repeated for subsequent steps 
5 until there is convergence 38 whereupon one exits this inner 
loop. 

Consider the outer loop. After exiting the inner loop 38 
the solution for {h f } is tested to see if it meets the 

10 performance goals 39* If not, a new set of metric weights is 
selected and the inner loop is initialized i=0. 40 and the 
inner loop is used to find the next solution for {h f }. This 
process continues until the performance goals are met or 
adequated approximated 39 whereupon the algorithm is exited 

15 with the final solution set {h f } 41. 

FIG. 5A,5B,. . .,,arc tho Matlab 5.0 code for the LS 
recursive eigenvalue solution algorithm used to design — fcke 
Wa velet Wavelet in FIG. 6 . Included are a listing of tho frcq teeftey 

20 domain design harmonics, timo response/ and data plotting 
software used in the optimization/ and an example of the scal -ift^ 
e# this mother Wavelet into a Wavelet with now scalin g 
-d ilation) / time translation/ frequency translation paramo ^e^e-r 
ee d the functions used in the code. The Matlab code us ee— aft 

25 iterative approach to find the optimal frequency d emajrft 
ei genvector solution to tho quadratic form for tho LS error whos e 
error matrix is the weighed sum of the error matrices for -4he 
efer opband/ passband/ dcadband metrics and the error matrices fog 
^r ko ISI and ACI metrics which arc functions of the cigcnvcct eg3-=- 

30 For each iteration tho ISI and ACI error matrices arc constr tteted 
with the frequency domain eigenvectors from the prcvioa g- 
it eration. The iteration finds the new frequency domain 
ei genvectors which arc then used to update the ISI and — A€hE 
metrics for the next iteration. This iteration is repeated until 

35 there is convergence. The Matlab code allows tho operator ^fee 
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update the metric woighto ao well ao the oGcnario paramctGra — fee- 
find the optimal choice for the performance goaJ r^-r 

Matlab code otarto with the oclcction of the dc& jr€fft 
5 paramctcro in 29 Stop 1, Stop 1.1 lioto the occnario paramet ers- 
wk ioh are digital oamploo per Wavelet ly oamplo intorval r 

L-^16 nominal length of Wavelet y in unito of M, N~ML4l Wavelet iy 
jre ngth which io rcprcocntcd by N' in oquationo (10) , fo~l i -s-^fche 
fto rmalieod channel opacing equal to the Wavelet \y sample ratc r 

10 fp-0.8864 io the paooband frequency interval; n f~16 is — fche 
fttt mbcr of dcoign frequency harmonico, n fft~1024 io the number o f 
digital oamplco uocd for calculation of the frequency opectrum 
centered at 0 frequency, obno^fib/^^ . 0 dB io thodeoign val ge-ef 
the energy per bit E b - ^to noioc power denoity -&o—tttio , and 

15 x imbal aoi~6.0 dB io the aooumcd channel to channel power 
Afflb alancc uocd in the calculation of the ACI crroro. Stop 1.2 
de rivco the ooftwarc paramotoro from the dooign parametoro in 
S- feop 1.1. Stop 1.3 lioto the optimisation paramctero which -age 
fefe c number of itcrationo n iteration - 10 / and the metric wc i-ghfeg- 

20 w (metric) from cquat jreft — frS^j and in FIG. 4 equal to w(pao e)-= 
alpha 1-1 . c - 2 for paooband/ w (otop) -alpha 2-0.80 for otop feaftdr- 
w-e^S I) -alpha 3-2. c 3 for ISI, w (ACI) -alpha 4-0.5 for ACI, an d 
frfrf dead) -alpha 5~0 for dcadband. The number of iterations — wee 
se lected to provide a convergent oolution for the weighted L & 

25 error mctrico w (metric) J (metric) in FIG. SK 38 Step 10 figu ge-f£f 
and for their oum J in figure (1). Weight valueo wore ooloct ed-^fce 
optimise the Wavelet \\t filter performance in FIG. 6 and 
figurcQ) in FIG. 5L 3S Step 11/ Wavelet v ripple/ ISI/ ACE 
e-i gnal to - noioo SNR power ratio loooeo in figure (4)/ and Wav ele-fe 

30 \\f time rooponoo in figure (5 -K- 

¥iGnr SB 30 Stop 2 arc the initialization calculationo ^jre^ 
%e- the otart of tho iteration loop. Stop 2.1 calculateo thc -eee— 
ei- ded Wavelet \|/ length parameter m-128. Stop 2.2 calculatoo t he 
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matrix tranof ormation bw matrix which maps the \|/ frequency domain 
d esign cigchvcctoro into the one - aided \\f time response. Step 2.3 
-ce fcro to the function pmn in FIG. 5P,5Q 4 1 in Stop 13.1 which 
computes tho LS error matrix which io the sampled data integral 
5 £[8c » Se / ] do over the band in equation (23) for the the passba ftd 
a** d otopband metrics J (paso) > J (stop) . The function pmn d 
ee mputco the corrcoponding integral over the doadband in equation 
(2 4 ) and io not lioted since it wao not used for FIG. 6 : Step 
2. 4 constructs a structural matrix c matrix uocd in the ISI and 
10 AG I calculations. Step 2.5 constructs the sample rate templates - 
-£e r the pasoband upper cdgc>> sample rate marker/ and pasoban d 
lower edge for frequency performance evaluation. 

FIG. 5C 31 S tep 3 uses the function pmn to compute the LS - 
15 error matrices passband and otopband in Step 3.1, 3.2. Dcadb aftd 
h& error matrix dcadband is set equal to a null matrix in S top 
3.3. Step 3. 4 computes the matrix p total equal to the weighted 
-&u m of those matrices. Stop 3.5 maps this LS error matrix 
p— total in the time domain into the LS error matrix pw t in th e 
20 frequency domain using the matrix tranof ormation bw matrix. - 

FIG . 5C,5D 32 S tep 4 initialises matrices uocd in the 
-i teration loop and starts the iteration loop for tho iterations 

25 i iteration^!/.,., 10 . Stop 4.1 calculates the eigenvalue and 
e igenvector which minimize the cost function J in equations (27) 
aft d FIG. 4 whose LS error matrix io pw t. Stop 4 .2 uses tho 
b w matrix to map the eigenvector in frequency into the Wave jre^- 
time response ly-hn and otorco the eigenvector ao the set of 

30 Wavelet \|/ design harmonics — ^=h w oig(Ic) for Ic-Q, 1 /...-, 15 where \\f k 
h as boon defined in equations (9) and in equations (11) upon 
•ge cogniaing that there arc n f-16-^L design harmonics k^k^O,!, . y 
. /15. There io no roaoon to use tho negative harmonics since the 
spectrum is real and symmetric about the 0 frequency. Step 4-^3 

40 



computoo the paooband peak - to - peak ripple ao xripple and fehe 
otopband attenuation ao xot ep-r- 



5E 33 S top 5 caloulatco the weighted paooband cr^e^ 
me tric w (paoo) J(paoo) -beta paoo/ weighted otopband error motri e 
w(otop) J(otop)"bota otop/ and weighted dcadband error m etrE4re 
w-f dead) J (dead) ^bcta dead in cquationo (27) using the, b ve e^ef 
-ge rmat for hn and the LS error matrioco paooband^ otopba ftdr 
dcadband r 

FIC. 5E,5F 34 Step 6 calculatco the ISI and ACI LS error 
Hta tricco/ their mctrieo/ and their error contributationo t e-^fcke 
looo. Step 6 .1 oaleulatoo the ISI LS error matrix w matrix 
aft d the oorrcoponding ISI metric J (ISI) ~crrM ioi with Matlab co de 
15 which implemcnto cquationo (25) / and calculatco the ISI rcoidu ark 
e^ ror orrV ioi that cauoco an ISI SNR looo. Step 6 .2 calculat es- 
the ACI LS error matrix w f matrix and the oorrooponding — 
me tric J (ACI)- crrM aci with Matlab code which implemcnto 
e^ uationo (2 6 ) , and calculatco the ACI rcoidual error crr ¥— aejr 
20 that cauoco a ACI S NR looo. 

5C 35 Step 7 calculatco the weighted LS error mctrieo 
•fo g ISI and ACI and the LS error matrix pw t for the ne xt: 
jrfc eration. Stop 7.1 calculatco the weighted I S I LS error metric 
25 w(ISI) J(I S I)-bcta ioi and the weighed ACI LS error mct gjre 
w (ACI) J (ACI) "beta aci in cquationo (27) and the oum of the 
we ighted LS error motrico J-errM ioi. Step 7,2 oaveo the weighte d 
L& error mctrieo and their oum J for each iteration. Step 7.3 
updatco the pw t matrix for the next iteration. 

30 

-F£G-7- 5C,6H 3 6 Step 8 calculatco the SNR loooco, otor ee 
t& ooo loooco for each iteration/ and complctco the iteration 
jre op. Stop 8 .1 computoo the SNR looo in dB due to paooband 
« pplo ao xlooo ripple; due to ISI ao xlooo ioi/ due to ACI aa 
35 xlooo aci, and the total SNR looo ao xlooo total. Step 8.2 oaves 
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thooc SNR looooo for each iteration and oomplctco the itcrat -aref* 
loop. 



F£G^— 5H,5I,5J,5K 37 Stop 9 dooumcnto tho Wavelet \y 
5 frequency dooign harmonico and timo rooponoc goncratcd by tho 
Ma tlab codo for tho Wavelet in FIG, 6 . Step 9.1 lioto t he 
Wavelet vy frequency domain harmonic dcoign coordinates 
j fa-f k)-hw cig(k) for harmonico k-Q, I,.../ 15. Only tho pooitivc 
h armonico are lioted oinoc the frequency opectrum io real and 
10 oymmotric about k-^0. Stop 9.2 lioto the Wavolot vy time rco pof*se 
\H* Q-hn(n) for n-Q,l,2,. . . ,128 with the coordinate frame 
ee ntorcd at tho Wavelet peak. Only the pooitivc digital oa mpjre^ 
arc lioted oincc the time rooponoc io real and oymmctric ab e^t 
n^O. 

15 

. FIG. 5K 3 8 Step 10 ploto the total weighted LS error metric 
3— in figure (1) and tho individual weighted LS error metrico in 
#3r gurc(2) ao functiono of the iteration number. After about 6-7 
jet crationo both the total and the individual, weighted LS crror - 
20 metrico arc oboerved to otabiliac indicating that the oolcctcd 
member of . itcrationo n itcratiori-10 io oufficient for - 
convergence . 

FIG. 5L 39 S tep 11 ploto the Wavclot vj/ frequency rooponoc 
25 in figure (3) which io the oamc ao the froquoncy rooponoc in FIG. 
•frr- ploto the SNR loooco due to paooband ripple^ ISI, ACI and the 
%-e- tal loop in figure ( 4 ), and ploto the time rooponoc in figure (5 )- 
for uoo in the optimisation of the LS metric wcighto in Step 1.3 
ao well ao the occnario paramotcro in Stop 1.1. 

30 

FIG. 5M 4 0 Step 12 documcnto tho Matlab code to reocale t -he 
ae ther Wavclot y in FIG. 6 which wao derived in the previouo 
Matlab code in Stcpo 1 - 11> for new ocalc (dilation) p> time 
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tranolation q, and frequency tranolation k parameters comparcd --%e 
the valuoo p-Q, q-Q/ k-Q for the mother Wavolct \y. An example ee^ 
e-f new parameter o documented in the code io p-2> > q~2 / k-3 and 
wjr th the aooumption that the digital oampling rate remains fixe d 
5 under the parameter change. A conotant digital oamplc rate kecpo 
the frequency band conotant under the parameter change. 

FIG. 4,5 4012.1(7) , (1 8 ) (2 8f 

New Wavelet scaled from mother Wavolct (2 8 ) 

10 

^ P ,g,H,L.i^ e^ ^= ^D . q ,k with paramctero M, L 

- — gr* ^ y (2"P n ■ qM) cxp ( j2rc le-g-^n /ML) 

There arc at lcaot two mcthodo of ocaling uocd to implement 
15 equation (28) : In method 1 porrcoponding to the otandard Wavc jre^- 
aooumption, the ouboampling or decimation of n by the factor 
•& tretchoo the time between contiguouo oamplco by the factor 2^p 
wh ale keeping conotant the number of oamplco M per Wavelet y 
ea mplc interval/ and therefore roducco the frequency band by th e 
20 factor 2 A p. For a uniform filter bank/ thio mcano the fil feeg 
fe cquency bando arc reduced by the factor 2 A p. In method 2, the 
o ampling io held conotant corrcoponding to a conotant frequency 
ba fid aooumption. The factor 2 A ( - p) io divided out in the Wavel et 
e xprcooion which io equivalent to the M being ocalcd to 2 A pM 
25 corrcoponding to increasing the number of oamplco per Wavelet 
interval by the factor 2 A p. For a uniform filter bank/ thio 
me ano method 2 corrcopondo to incroaoing the maximum number o# 
eh annclo by the factor 2 A p over the oamo frequency band. 

30 

Caoo 1 Scaled (dilated) oampling 2~ * *n with fixed M 
— ouboamplco (decimatco) n by the factor 2 P 

■ Rcprcocnt n ao 
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n - n 0 i 2 P n p 

n 0 - 0,1, . . . , 2 *-± 

n p - 0, 1/ 1, i/ 2, . . . 

■ Wc obocrvc 2~^n - n p 

5 * Thio cnablco the now Wavelet to be written 

vi/ D ,Q,k ^ -Z^- Vin p qM) oxp(j2ic k n p /ML) 

Cqoo 2 Sampling n romaino fixed and M io 
repealed by the factor 2 P 
10 ■ New M io M new 

M now ^2^14 

■ Thio cnablco the now Wavelet to bo writt en: 

y/ 1 & rqHe -^-2-~ p ^ \\f (n - qM new) cxp(j2rc kn/M -ftewirf 

15 

wh erein tho new Wavelet equation from (7) , (1 8 ) io rc - writto Br 
uoing the frequency tranolation parameter k^-r and tho expli ejrfc 
frequency tranolation cxprcooion oxp ( j2rc k 2 - ^n /ML) which -is- 
20 recognized ao the diocrotc Fourier tranoform kornal. 

In Caoe 1 corrcoponding to the otandard Wavelet oooumptiony 
^h:C ouboampling or decimation of n by the factor 2 A p otrctchc e- 

25 the time between contiguouo oamploo by the factor 2 A p wh irle 
■k eeping conotant the number of oamplco M per Wavelet \y oamplc 
•ifi- torval/ and therefore roducoo tho frequency band by the factor 
2-^p . Tho now n can be rcprcoontcd ao tho oum of tho two field s 
fl-G- r 0,1,. . , (2 A p 1) and n p - 0, 1 / 2, . . . Scal -jrtt^ 

30 or oquivalcntly ouboampling rooulto in the new n equal to n -new 
— 2 A - p n - n p corrcoponding to a otrctching of the oampl e 
jbft tcrval by 2 A p. The new Wavelet length N new rcmaino conotant 
N now~N uoing tho n p oampling. 
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In Caoc 2 which io the communicationo application being 
e onoidcrcd, the oampling io held conotant corrcoponding to a - 
conotant frequency band aoaumption. The factor io divi ded 

e* *t in the Wavelet cxprcooion which io equivalent to the M being 
5 ocalod to M now-2 A p M corrcoponding to incrcaoing the number of 
s amplco per Wavelet interval by the factor 2 A p. Thio mcano Caoc 
2 corrcopondo to incrcaoing the maximum number of channolo by ^he 
-f actor 2 A p over the oamc frequency band. The new Wavelet length 
jre- equal to N ncw~M nowL + 1 -1025 oince M ncw~2 A pM^&4-^ 

10 

FIG. 5M 4 0 S top 12.2 dcrivco the new matrix tranof ormation 
b w matrix new which mapo the \y frequency _ domain dcoign 
ei genvectoro into the onc - oided i|/ time rcoponoo. 

15 FIG. 5M,5N 40 Stop 12.3 uoco the matrix trano format i -eft 

bw matrix new to map the deoign eigenvector hw cig into the new 
Wa velet in three otcpo. In the firot otcp> the mother Wavelet io 
^anolatod into the new baocband Wavelet -* p P , q - o,k - o ~ hn 0 which io 
£h-c mother Wavelet ocalod by p and with no time and frequency 

20 tranolationo q-lc-Q. The oocond otep tranolatc thio baocband 
Wavelet by ,qM new to generate -*j * p , q ,ic - o " hn 1 . The laot otep 
tranolatco thio Wavelet M/ p , qv k - o - hn 1 in frequency to give the -new 
Wavelet — *f f p , q ,k =• hn new. 

25 FIG. 5N 4 0 Stop 12. 4 ploto the Wavelet time rooponoeo for 



the — mother 


Wavelet — 


4m — a**d — fcke — aew — Wavelet — *p p a , k = — hn now 


— if* 


figure (6) . 






for 




Stop 12.5 


ploto the Wavelet frequency rcoponooo 




the mother 






in 




Wavelet — 


4«i — aftd — the — now Wavelet — *p p , q ,) C = — hn new 




figure (7) — vorouo the normalized frequency/ (ly-hn oamplc rate) 


and 














in figure — 


-f8-) — vorouo 


the — normalized — frequency/ (i|/ p a k = — 


now 



oamplc rato) . 
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FIG. 4 1 S top 13 dooumonto the functiono uocci in tho Matlab 
code. Stop 13.1 io the code for the function pmn which computet? 
fe-h e matrix for the J (BAND) in cquationo (23). FIG. 4 1 Stop 13,2 
j^r oto the code for the function frcq rop which computco the 
5 Fourier tranoform of tho input' hn vcrouo tho normali eed 
-g^ oquoncy/ (Wavelet y-hn oamplc rate). 

Conoider the inner ncotcd loop 2 8 . The recuroivo L& 
se lution otarto with the initial otep irQ 29 which begins with 

10 tho ooloction of tho input dcoign paramctoro and wcighto and th e 
se lection of an initial oct of valuco for the (h^) 30 . Next the 
ba nd mctrico are calculated 31. S ince thio io the initial otep 
4r =0 32 the coot function J io rcotrictcd to the linear band 
me trics in order to fine the linear approximation to [h^ f 33 t e 

15 initialize the non-linear oolution otarting with otep i~l —34 
w herein tho highly non linear interference mctrico arc included 
•ana- J. Following tho oignal flow> for otep i~l the band metrics - 
33r and interference mctrico 35 are calculated/ weighted^ and 
ex amined to form tho coot function J 3 6 . An iterative oolution 

20 algorithm 3? findo the boot approximation in otep i-1 to the 
-fkjj which minimizco J uoing the approximation to [h#} from tho 
p-p eviouo otop i-Q to linoarigc tho ocarch algorithm cocfficicnto 
£er the eigenvalue and LS otep i~l iteration. If there io n & 
ee nvorgence to tho correct oolution for — — thio rccuroiv e 

25 ocquoncc of calculationo io repeated for otep i~2 3 4 . Thio 
recur oivc oolution tochniquo io repeated for oubooquont otopo 
ttft til there io convergence 38 whereupon one cxito thio inner 
loop. 

Conoider the outer loop. After exiting tho inner loop 3 8 
30 the solution for — frh *) io tooted to ooo if it moeto fe he 
performance goalo 39. If not> a now oct of metric wcighto i& 
•se lected and tho inner loop io initialized i^O 4 0 and th e 
jHft ncr loop io uocd to find the next oolution for [h^ H- Thi& 

proccoo' continuoo until the performance goalo are met or 

■ J 



46 



ad cquatcd approximated 39 whereupon the algorithm io exited 
with the final oolution oet [hd — 4 1, 

Applications of this new invention to both communications 
5 and radar will be given using these example algorithms and other 
algorithms supported by this invention. These new waveforms are 
considered for the applications: l)to replace the square-root 
raised cosine waveform (sq-root rc) which is extensively used for 
the third generation (3G) CDMA communications, 2) to replace the 
10 Gaussian minimum shift keying (GMSK) waveform for constant 
amplitude bandwidth efficient (BEM) applications, and 3) as a 
candidate waveform for synthetic aperture radar (SAR) and real 
aperture radar (RAR) applications. 

15 CDMA communications application for the? current and the new 

3G CDMA considers a waveform designed with this new invention as 
a possible replacement for the sq-rt rc waveform with bandwidth 
expansion parameter a = 0.22 to a = 0.40. This notation means that 
for a =0.22 the spectral efficiency is (symbol rate/bandwidth) = 

20 1/1 + a =1/1.22 = 0.82 = 82%. A basic advantage of the waveform 
is the potential for a symbol rate increase within the same 
bandwidth with an increase in the spectral efficiency to «100% 
, depending on the application and operational constraints. The dc 
power spectral density or power spectrum (PSD) of the waveform is 

25 compared to the PSD for the sq-rt r-c in FIG. €8. Plotted are the 
measured PSD in dB 42 versus the frequency offset from dc 
expressed in units of the symbol rate 43 Plotted against the 
normalized frequency offset are the dc PSD for the new waveform 
44, the sq-rt r-c with <x=0.22 45, and the sq-rt r-c with 

30 a=0.40 46. It is observed that the PSD for the new waveform 
rolls off faster that that for the sq-rt r-c which means that the 
new mr waveform will support an increased symbol rate for a given 
available frequency band while satisfying the inherent 
requirements for low ISI and MAI (multiple access interference) . 
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Constant amplitude BEM application of the new waveform 
indicates that it is a viable candidate for replacing the current 
preferred modulation waveform which is the GMSK. The GMSK finds 
5 applications for transmitters which operate their HPA(s) 
amplifiers in a saturation mode in order to maximize, their 
radiated power from the HPA(s), and which require a BEM PSD to 
avoid excessive spreading of the transmitted power. Simulation 
data for the dc PSD is plotted in FIG- 3—9 for the new waveform 

10 BEM and the GMSK. Plotted are the measured dc PSD in dB 47 
versus the frequency offset from dc expressed in units of the bit 
rate 48. Plotted against the normalized frequency offset are 
the dc PSD for the new waveform BEM 49 and the GMSK 50, for a 
length parameter L=10 where L is the length of the phase pulse in 

15 terms of the phase pulse repetition rate. The significance of 
this example data is that the new waveform has the potential to 
be designed to offer a PSD which is less spread out than the 
current GMSK, and therefore an improved BEM waveform. 

20 Real aperture Radar — radar ( RAR) and synthetic aperture 

radar JSAR_)_ applications of the waveform indicates that it is a 
viable candidate to replace the current chirp waveforms for 
wideband signal transmission, when combined with pseudo-random 
phase codes. Results of the simulation for the new waveform and 

25 an unweighted frequency chirp waveform are given in FIG. -8-10. 
Plotted are the ambiguity function for the new waveform 51 and 
the unweighted frequency chirp waveform 52. The dc 2- 
dimensional radar ambiguity function 53 is plotted as a function 
of the frequency offset in units of fT p and the time offset in 

30 units of. t/T c where T p is the phase-coded radar pulse length or 
length of the phase code and T c is the phase code chip length. 
The chip length is identical to the waveform repetition interval 
T s so that T c = T s . It is observed that the new waveform has the 
potential for significant improvements in the ambiguity function 

35 and by implication in the performance. 
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Preferred embodiments in the previous description is — are 
provided to enable any person skilled in the art to make or use 
the present invention. The various modifications to these 
5 embodiments will be readily apparent to those skilled in the art, 
and the generic principles defined herein may be applied to other 
embodiments without the use of the inventive faculty. Thus, the 
present invention is not intended to be limited to the 
embodiments shown herein fe^band— is not to be accorded the wider 
10 scope consistent with the principles and novel features disclosed 
herein. 
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PIG. 1 Wavelet Tiling of an N- Point Di/ital t-f Space 
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FIG. 2 Wavelet Iterated Filter Bank for Tiding t-f 
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FIG. 3 PSD Requirements for Communications 




FIG. 4 LS Metrics and Cost Function 
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FIG. 5 LS Recursive Solution Algorithm 



Start 



29 



i=0- 



40 



30 



Select Design Parameters 



31 



35 
36 
37 

38 



Calculate Band Metrics 



:;;-...:, ; ^32<i= 



39 



Solution 




v Reset i=0 ' J:-- - 
• Select new {w"(me tries).} 



new {h f } v 7 "' 

* i=i+l . . '.. * 



ffi&ft 



Find {h f } ,at step'- i==0 
to minimize- cost - J " 



jbnvergenc^fto solution 



: Inner, 'Loop::: 
Recur si ve LS > 

.... ... * ■ - ^ ■ •:*' •'• „■ tv\.'-r . ■ „ ■ " ■:, ... 

St' J? v . -v.^^ 

ii : \:'" ; <: 7 *. ^XV^^v * * W'^f^ 



goals -?■■ 



Solution set 



Outer 
Loop 27 



55 



FIG. 6 PSD for New Waveform and Square-Root Raised-Cosine 




FIG. 7 PSD for New Waveform BEM and GMSK 
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